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ABSTRACT 


The  Center  of  Excellence  in  Hypersonics  Research  was  established  to  support  future  Department 
of  Defense  and  NASA  hypersonic  vehicle  development  programs.  The  Center  engaged  leading 
experts  in  experimental  and  computational  analysis  of  hypersonic  flows  to  provide  research, 
education,  and  technology  development  in  hypersonics.  The  Center  combined  the  world’s 
premiere  hypervelocity  testing  facility  at  the  Calspan  -  University  at  Buffalo  Research  Center 
with  state-of-the-art  computational  simulation  methods  at  the  University  of  Minnesota  to  rapidly 
advance  the  field  of  hypersonics.  The  Center  supported  fundamental  research  in  hypersonics, 
including  the  development  of  high-quality  experimental  databases  for  code  testing  and  validation 
of  prediction  tools.  Many  of  these  codes  have  been  transferred  to  government  and  industry  to 
enable  optimized  hypersonic  vehicle  designs.  The  Center  also  supported  a  rigorous  graduate 
program  that  transferred  technical  capabilities  from  the  previous  generation  of  hypersonics 
experts  to  new  students  in  the  field.  This  helped  preserve  the  nation’s  dwindling  talent  pool  and 
provided  a  source  of  scientists  and  engineers  for  future  research  and  development  programs.  The 
Center  of  Excellence  in  Hypersonics  Research  was  a  national  focal  point  for  U.S.  research  and 
training  in  hypersonics. 
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INTRODUCTION 


The  major  programs  currently  being  undertaken  by  the  Department  of  Defense  and  NASA  to 
develop  new  advanced  hypersonic  vehicles  and  space  access  systems  will  require  significant 
advances  in  the  design  methods  and  ground  testing  techniques  to  ensure  successful  flight  testing. 
In  the  period  from  1950  to  1970,  the  U.S.  led  the  world  in  the  development  of  space  access 
systems  and  hypersonic  vehicle  design.  However,  the  past  decade  has  seen  the  loss  and 
retirement  of  many  of  the  talented  personnel  trained  in  this  era  and  a  dearth  in  research  to  resolve 
key  problems  associated  with  aerothennal  problems  arising  in  hypersonic  vehicle  design. 
Exacerbating  this  problem  has  been  the  lack  of  funding  to  train  new  scientists  in  this  area  and  the 
development  of  numerical  and  experimental  techniques  to  provide  accurate  models  of  the 
aerothermodynamic  processes  used  in  the  methods  to  design  and  analyze  hypersonic  vehicles 
employing  revolutionary  designs.  Such  is  the  status  of  this  area  that  there  is  significant  question 
as  to  whether  currently  America  has  the  engineering  talent  to  rapidly  replicate  the  activities  in 
the  1960s  which  placed  American  astronauts  on  the  moon.  These  observations  motivated  the 
need  for  a  Center  of  Excellence  in  the  area  of  hypersonic  flows  to  advance  the  state  of  the  art  in 
the  understanding  of  hypersonic  flows  over  vehicles  and  to  train  the  scientists  and  engineers  who 
will  design  these  future  hypersonic  vehicles. 

Under  AFOSR  funding,  we  established  a  partnership  between  the  University  of  Minnesota  and 
the  Calspan  -  University  at  Buffalo  Research  Center  (CUBRC)  to  develop  a  Center  of 
Excellence  in  Hypersonic  Research.  CUBRC  has  the  world’s  premiere  hypersonic/hypervelocity 
testing  facilities  and  unique  hypervelocity  instrumentation  developed  over  40  years  of  wind 
tunnel  testing  at  hypersonic  flight  conditions.  The  University  of  Minnesota  has  state-of-the-art 
capabilities  to  simulate  hypersonic  flows  including  laminar  complex  regions  of  shock/boundary 
layer  interaction  and  flow  chemistry,  transition  to  turbulence,  finite-rate  air  and  combustion 
chemistry,  turbulent  mixing  in  compressible  flows,  and  shock  wave/turbulent  boundary  layer 
interactions.  The  research  supported  by  AFOSR  through  the  Center  has  greatly  advanced  the 
ability  to  accurately  test  and  simulate  hypersonic  flows. 

The  Center  trained  many  young  engineers  and  scientists  to  conduct  experimental  research  in 
state-of-the-art  hypersonic  test  facilities  at  CUBRC,  and  perfonn  numerical  simulation  and 
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analysis  with  the  most  advanced  numerical  codes  and  modeling  techniques  at  the  University  of 
Minnesota.  A  major  research  focus  was  a  program  of  combined  numerical  and  experimental 
analysis  to  investigate  key  phenomena  associated  with  the  performance  of  re-entry  vehicles, 
airbreathing  propulsion  systems,  and  hypersonic  interceptors.  Employing  its  unique  combination 
of  large-scale  test  facilities  and  efficient  computational  simulation  methods  the  Center  stimulated 
advances  in  the  aerothermodynamics  of  hypervelocity  flows,  leading  to  improved  methods  to 
design,  analyze  and  evaluate  the  perfonnance  of  new  hypersonic  vehicles  and  airbreathing 
propulsion  systems. 
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OBJECTIVES  AND  APPROACH 


The  objective  of  the  Center  of  Excellence  in  Hypersonics  was  to  develop  a  strong  research  and 
training  program  in  hypersonics  to  support  the  Nation’s  future  hypersonic  vehicle  development. 
The  Center  drew  on  the  extensive  experience  of  the  principal  investigators  in  hypersonic  ground 
testing,  numerical  simulation,  graduate  student  advising,  and  classroom  teaching.  The  Center 
was  a  national  focal  point  for  academic  research  and  training  in  the  U.S. 

The  Center  of  Excellence  in  Hypersonics  Research  combined  fundamental  research,  education, 
and  technology  development  to  support  the  nation’s  future  hypersonic  vehicle  development.  The 
primary  mission  of  the  Center  was  to  perform  research  on  critical  problems  in  the  field,  such  as 
real  gas  effects,  shock  wave  -  turbulent  boundary  layer  interactions,  transition  to  turbulence,  and 
scramjet  testing  and  simulation.  The  research  program  was  used  to  drive  the  academic  mission  of 
the  Center,  which  was  to  train  a  new  generation  of  hypersonics  experts.  This  was  accomplished 
through  relevant  course  work,  hands-on  computational  method  development,  design  and  testing 
of  experiments,  and  computational  analysis  of  experimental  results. 

KEY  ACCOMPLISHMENTS 

The  AFOSR  Center  of  Excellence  in  Hypersonics  Research  was  a  large-scale  multi-year  effort 
that  resulted  in  a  very  large  number  of  publications;  these  are  listed  below.  It  is  not  possible  to 
provide  detailed  infonnation  about  each  sub-topic  of  the  research  accomplished  under  this 
program.  Rather,  in  this  section  we  provide  key  highlights  of  the  research  and  refer  the  reader  to 
the  publications  for  additional  details.  We  organize  the  work  into  several  major  topics: 

1.  Development  of  advanced  numerical  methods  and  codes  for  high-speed  and  hypersonic 
flow  simulations. 

2.  Development  of  two-dimensional,  axisymmetric  and  fully  three-dimensional  boundary 
layer  stability  analysis  tools. 

3.  Detailed  experimental  measurements  to  provide  validation  data  for  hypersonic  flow 
simulations. 

4.  Hybrid  RANS-LES  simulations  and  direct  numerical  simulations  of  transitional  and 
turbulent  hypersonic  flows. 
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5.  Advanced  particle-based  simulation  methods  (Direct  Simulation  Monte  Carlo  and 
molecular  dynamics)  for  hypersonic  rarefied  flows 

Each  of  these  topics  has  produced  many  papers  and  several  Ph.D.  theses  over  the  course  of  the 
Center  program.  In  many  ways,  this  work  leads  the  field  and  it  would  not  have  been  possible 
without  the  AFOSR  support  of  the  Center.  Much  of  this  work  has  been  transitioned  to  Air  Force 
Research  Faboratories  and  to  other  national  labs  and  industry. 

Advanced  Numerical  Methods  for  Hypersonic  Flow  Simulations 

The  standard  numerical  simulation  approach  for  high-speed  shock-dominated  flows  has  been 
based  on  upwind-biased  Riemann  solvers  since  the  early  days  of  computational  aerodynamics. 
These  numerical  flux  functions  are  very  effective  at  providing  stable,  robust  simulations  of 
steady-state  hypersonic  flows  including  the  effects  of  chemical  reactions.  However,  these 
methods  stabilize  the  flow  by  adding  large  levels  of  numerical  dissipation  -  simply  put,  the 
upwind  flux  functions  smear  out  discontinuities  and  cause  small-scale  unsteady  motion  (kinetic 
energy)  to  be  converted  to  thennal  energy.  Thus,  when  these  methods  are  used  to  simulate 
unsteady  flows,  they  artificially  dissipate  the  small-scale  (both  in  length  and  time)  unsteadiness. 
To  accurately  resolve  an  unsteady  flow  with  these  methods  would  require  unpractically  large 
grids,  especially  as  the  Reynolds  number  approaches  realistic  flight  conditions.  Therefore,  two 
new  classes  of  methods  were  developed  under  Center  support. 

Two  approaches  were  taken.  In  the  Candler  group,  a  new  numerical  flux  function  was  developed 
that  conserves  a  secondary  quantity,  in  this  case  kinetic  energy  [63].  For  a  compressible  flow, 
kinetic  energy  is  not  actually  conserved,  so  for  compressible  flows  this  approach  is  consistent 
with  kinetic  energy  conservation.  Namely,  the  resulting  discrete  flux  of  kinetic  energy  is 
consistent  with  the  kinetic  energy  conservation  equation.  This  approach  was  shown  to 
significantly  reduce  the  numerical  dissipation  for  standard  numerical  test  cases.  For  example,  the 
method  produces  the  expected  behavior  as  a  function  of  Reynolds  number  for  the  decay  of 
weakly  compressible  and  highly  compressible  isotropic  turbulence.  With  the  careful  addition  of 
dissipation  to  control  behavior  near  strong  shock  waves  and  density  gradients,  the  kinetic  energy 
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consistent  method  has  been  applied  to  a  wide  array  of  practical  problems.  These  will  be 
discussed  in  more  detail  below. 

Recently,  the  kinetic  energy  consistent  fluxes  have  been  extended  to  fourth  and  sixth  order 
accuracy  for  certain  classes  of  grids.  The  theoretical  order  of  accuracy  can  be  achieved  on 
weakly  stretched  hexahedral  grids,  as  shown  by  numerical  test  cases.  For  practical  problems  that 
involve  non-ideal  grids,  the  higher-order  formulations  have  been  shown  to  produce  clearly 
superior  solutions  than  the  lower  order  method.  Thus,  even  though  ideal  performance  is  not 
always  obtained,  the  higher-order  methods  are  very  useful  for  realistic  problems;  in  many  ways, 
they  enable  new  classes  of  simulations  of  turbulent  high-speed  flows. 

Work  was  also  conducted  to  improve  the  time  integration  aspect  of  numerical  simulations.  A 
novel  second-order  modified  Crank-Nicholson  approach  was  developed  and  tested.  For  certain 
problems  (e.g.  jets  in  supersonic  crossflow),  the  method  extends  the  accuracy  of  the  simulations 
allowing  larger  time  steps  to  be  taken  for  a  given  level  of  accuracy. 

In  the  Mahesh  group,  a  new  approach  that  rescales  the  governing  equations  so  that  they  become 
inherently  less  stiff  (difficult  to  solve  numerically)  at  low  Mach  numbers,  yet  still  producing 
accurate  results  at  high  Mach  number  [17].  This  approach  is  documented  in  the  cited  Journal  of 
Computational  Physics  articles  and  has  been  applied  to  numerous  compressible  flow  problems. 

Both  approaches  have  been  implemented  in  large-scale,  highly-parallel  computational  fluid 
dynamics  codes.  The  US3D  code  [5,  11]  has  been  transitioned  to  several  AFRL  sites  including 
Wright-Patterson  AFB  (Air  Vehicles  Directorate)  and  Kirtland  AFB  (Space  Vehicles 
Directorate).  This  code  includes  many  key  features  for  practical  flow  simulations:  general 
thermo-chemical  model  for  high-enthalpy  nonequilibrium  flows,  laminar  and  RANS  transport 
models,  capability  of  rapid  convergence  to  steady-state  using  implicit  methods  and  upwind 
fluxes  or  time-accurate  mode  with  low-dissipation  kinetic  energy  fluxes  and  higher-order  time 
accurate  integration  methods,  the  ability  to  handle  hybrid  unstructured  grids  composed  of  a  mix 
of  tetrahedral,  pyramids,  prisms,  and  hexahedral  grid  elements.  (Though  hexahedral  grids  are 
generally  preferred  for  accuracy  reasons.)  A  critical  aspect  of  the  work  involved  designing  the 
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code  to  be  scalable  to  large  numbers  of  processing  elements  (CPUs  or  CPU  cores)  so  that  very 
large  simulations  can  be  performed.  This  required  careful  attention  be  paid  to  input/output 
processes  which  can  be  a  serious  limitation  for  large  grids;  the  code  reads  and  writes  in  parallel 
and  does  dynamic  partitioning  to  make  it  efficient  and  flexible.  It  has  been  run  successfully  on 
4000  processing  elements  on  a  500  million  element  grid.  There  is  no  inherent  reason  that  it 
should  not  scale  beyond  this  problem  size. 

Development  of  Tools  for  Hypersonic  Boundary  Layer  Stability  Analysis  and  Prediction 

A  key  accomplishment  from  the  Center  was  the  development  of  the  STABL  suite  of  codes  for 
hypersonic  boundary  layer  stability  analysis  (STABL  =  Stability  and  Transition  Analysis  for 
hypersonic  Boundary  Layers)  [12,  24,  60].  The  primary  component  of  the  STABL  code  is  a 
parabolized  stability  equations  (PSE)  solver  for  high-enthalpy  flows;  this  code  can  be  run  in 
parallel  to  reduce  simulation  times  and  it  can  be  driven  with  an  intuitive  graphical  user  interface 
so  that  non-experts  can  run  analyses  and  interpret  the  data.  It  includes  a  full  thermo-chemical 
model  for  nonequilibrium  flows  so  that  the  effects  of  finite-rate  chemical  reactions  and  internal 
energy  relaxation  can  be  included  in  the  stability  analysis.  In  addition,  STABL  includes  a 
parallelized  computational  fluid  dynamics  code  that  is  set  up  for  rapid  mean  flow  solutions  to 
support  the  stability  analyses.  However,  STABL  is  limited  to  two-dimensional  and  axisymmetric 
problems;  though  it  has  been  run  successfully  on  symmetry  planes  of  more  complex  three- 
dimensional  hypersonic  flows.  STABL  is  in  wide  use  by  US  national  labs,  academia,  and 
industry;  it  is  the  de  facto  standard  for  stability  analyses  of  hypersonic  flows.  It  was  used  to 
support  the  design  of  the  AFRL  HIFiRE  Flight  Experiments  (flights  1  and  5)  [34,  46-48,  50,  69], 
to  analyze  the  X-43  post-flight  data,  to  provide  predictions  of  transition  to  turbulence  on  the 
DARPA  HTV-2  re-entry  vehicle  [56],  and  it  is  has  been  and  is  currently  being  used  to  correlate 
Sandia  National  Laboratories  flight  data,  among  other  programs.  An  example  of  the  utility  of  the 
STABL  code  is  shown  in  Fig.  1,  which  plots  the  predicted  transition  location  versus  the 
measured  transition  location  for  two  prediction  methods  -  the  conventional  correlation  of 
Reynolds  number  based  on  momentum  thickness  over  boundary  layer  edge  Mach  number  and 
the  present  mechanism-based  approach  the  uses  stability  theory  to  compute  the  growth  of 
disturbances  and  correlates  their  integrated  growth  to  wind  tunnel  data  (the  <?A  method).  This 
result  is  for  the  HIFiRE- 1  blunt  cone  configuration  in  the  CUBRC  and  NASA  Langley  wind 
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tunnels.  Note  that  the  traditional  approach  has  very  limited  predictive  capability,  while  the 
mechanism-based  approach  correlates  the  diverse  data  very  well. 


(a) 


(b) 


Figure  1.  Predicted  vs.  measured  transition  location  on  the  HIFiRE-1  blunt  cone  using  two 
approaches:  (a)  the  mechanism-based  STABL  calculation  and  the  N-factor  approach  with  N  = 
5.5;  (b)  the  Ree  / Me=  150  empirical  approach. 
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A  three-dimensional  version  of  STABL  was  developed  under  Center  support.  This  code  required 
a  complete  reformulation  and  rewrite  of  the  stability  solver.  It  is  now  being  extended  to  include  a 
graphical  user  interface  and  features  that  will  make  is  useful  to  the  user  community. 

Experimental  Data  for  Code  Validation 

An  important  and  extensive  component  of  the  Center  was  collaborative  computational  and 
experimental  research  to  develop  high-quality  validation  data  for  high-enthalpy  flow  simulation 
codes  [1,  3,  5,  7,  8,  10,  32,  39,  83,  84],  This  work  largely  focused  on  the  double-cone  and 
hollow-cylinder- flare  configurations  that  have  been  the  subject  of  code  validation  studies.  In 
addition,  the  test-section  conditions  in  the  CUBRC  LENS-I  reflected  shock  tunnel  facility  were 
characterized  at  a  new  level  of  detail.  Laser  diagnostics  systems  were  used  to  measure  nitric 
oxide  (NO)  levels  to  be  measured  in  the  LENS  facilities  test-sections  [9,  22,  40].  This  work  was 
tightly  coupled  with  numerical  simulations  to  maximize  the  benefit  of  the  diagnostics  work. 
Critical  new  results  were  obtained  for  mixtures  of  argon  and  oxygen  at  a  range  of  enthalpies  for 
the  double-cone  flow  (discussed  in  more  detail  below).  These  measurements  show  that  the  lack 
of  agreement  between  simulation  and  experiment  at  high  enthalpy  is  due  to  the  dissociation  of 
oxygen  in  the  facility  reservoir;  present  thermo-chemical  models  for  the  expansion  of  this 
reacted  gas  through  the  facility  nozzle  are  inadequate  and  must  be  missing  some  key  physical 
process  or  processes.  In  particular,  it  is  possible  that  the  low-lying  electronically  excited  states  of 
oxygen  may  play  a  role  in  the  recombination  process;  work  is  continuing  in  the  attempt  to  better 
model  the  recombination  process  [80,  85].  In  the  following,  we  highlight  some  results  from  the 
recent  double-cone  experiments  and  related  simulations. 

High  Enthalpy  Measurements  on  25°/55°  Double  Cone 

In  response  to  observed  limitations  in  correctly  computing  the  laminar  separated  flowfield  over  a 
25755°  double  cone  geometry  (shown  in  Fig.  2),  a  number  of  new  experiments  were  perfonned 
targeting  two  areas  of  particular  difficulty  regarding  the  ability  of  simulations  to  reproduce 
experimental  measurements.  The  first  set  of  experiments  targeted  a  computational  unsteadiness 
observed  at  selected  high  enthalpy  nitrogen  conditions.  The  second  area  responded  to  limitations 
in  the  modeling  of  high  enthalpy  air  chemistry  by  generating  additional  data  using  mixtures  of 
oxygen  and  argon  as  a  test  gas,  thereby  eliminating  the  Zeldovich  exchange  processes  (involving 
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nitric  oxide)  from  the  physics  of  the  flow  field.  The  test  conditions  for  these  two  series  of  tests 
are  listed  in  Table  1. 


Figure  2.  Flowfield  over  25°/55°  double  cone 

During  the  course  of  a  recent  RTO  validation  study,  Run  40,  a  double  cone  test  case  run  at  high 
Reynolds  numbers  under  high  temperature  conditions  with  a  nitrogen  flow,  was  selected  for  code 
validation  studies.  In  retrospect,  this  choice  proved  unfortunate  because  while  the  similar 
measurements  in  air  were  computed  successfully,  and  a  stable  flow  was  predicted,  in  our 
program  no  such  computations  were  attempted  for  the  equivalent  run  conducted  in  nitrogen  - 
Run  40.  The  heat  transfer  and  pressure  measurements  made  over  the  double  cone  configuration 
in  Run  40  demonstrated  that  a  steady  separated  region  was  formed  within  4  ms  of  the  flow 
initiation  over  the  model  (Figure  3(a))  and  remained  steady  for  another  4  ms  as  illustrated  in 
Figure  3(b).  Figure  3(b)  illustrates  that  a  large  separated  region  is  fonned  at  the  cone/flare 
junction  with  a  well-defined  heating  peak  in  the  reattachment  region  followed  by  a  large  region 
well-defined  attached  flow  at  the  trailing  edge  of  the  flare.  There  are  no  indication  from  the 
surface  heat  transfer  and  pressure  measurements  or  from  the  high-speed  video  that  the  separated 
region  is  increasing  in  size  beyond  its  proportions  at  12  ms.  However,  for  this  test  condition,  a 
number  of  computationalists  have  obtained  solutions  where  the  separated  region  increases  in  size 
with  time  despite  the  fact  that  experiment  showed  no  such  flow  instability. 
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Figure  3.  Temporal  variation  of  the  heating  distribution  during  separated  flow  establishment  at 

the  double  cone  junction  for  Run  40 
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The  detailed  investigation  by  Nompelis,  et  al.  [84]  revealed  that,  for  his  computational  method, 
the  boundaries  between  stable  and  unstable  solutions  could  be  plotted  in  terms  of  the 
characteristic  Reynolds  number  and  the  freestream  total  enthalpy,  as  shown  in  Figure  4,  which 
indicate  that  decreasing  the  Reynolds  number  or  increasing  freestream  enthalpy  results  in  a 
predicted  stability  of  the  separated  region  across  a  well-defined  bounding  curve. 


Figure  4.  Nompelis’s  calculations  defining  the  boundaries  for  stable  and  unstable  N-S  solutions 

Based  on  these  calculations  we  selected  additional  high-enthalpy  test  cases  with  a  nitrogen 
freestream  flow.  From  these  runs  we  selected  (Run  80)  which  was  run  at  a  total  enthalpy  of  5 
MJ/kg  and  a  unit  Reynolds  number  half  that  of  Run  40  to  be  the  new  code  validation  test  case. 
Measurements  of  the  development  of  the  separated  region  for  this  test  case  is  illustrated  with  the 
heat  transfer  measurements  made  in  the  interaction  region  over  the  fore-  cone  in  Figure  5(a).  For 
this  condition  the  flow  is  established  in  the  test  section  in  approximately  1.25  ms  and  the 
separated  region  over  the  model  reaches  a  stable  size  in  approximately  4  ms  after  that.  As 
illustrated  in  Figure  5(b),  the  separated  region  remains  stable  in  size  for  the  following  5  ms 
during  the  remaining  steady  test  flow  through  the  tunnel. 
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Figure  5.  Temporal  variation  of  the  heating  distribution  during  separated  flow  establishment  at 
the  double  cone  junction  for  Run  80 
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The  predictions  made  by  Nompelis,  (which  are  similar  to  those  obtained  by  MacLean  at  CUBRC 
with  the  DPLR  code),  reach  steady  levels  in  approximately  the  same  flow  times  as  measured  by 
experiment  and  do  not  indicate  that  the  separated  region  grows  further  with  increasing  time. 
Comparisons  between  the  measurements  and  Nompelis  pre-test  predictions  of  heat  transfer  and 
pressure  through  the  interaction  regions  are  shown  in  Figure  6.  Flere  we  see  excellent  agreement 
between  prediction  and  measurements.  Clearly  this  test  case  (Run  80)  should  provide  a  more 
viable  code  validation  data  set  than  Run  40  for  CFD  validation  and  development. 


(a) 
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(b) 

Figure  6.  Steady  state  distribution  of  pressure  and  heat  transfer  over  the  double  cone 

A  serious  limitation  in  the  capability  to  predict  high  enthalpy  flow  over  the  double  cone  using 
reacting  air  as  a  test  medium  was  shown  by  Center  researchers.  In  response  the  these  findings,  a 
series  of  new  experiments  have  been  performed  to  eliminate  the  uncertainty  associated  with 
modeling  the  formation  of  nitric  oxide  where  oxygen  and  oxygen/argon  mixtures  were  used  as  a 
test  gas.  These  tests  used  pure  oxygen  at  several  enthalpy  levels  and  mixtures  of  various 
percentages  of  oxygen  diluted  with  argon  to  isolate  the  kinetics  of  oxygen  dissociation  alone. 

The  comparison  between  theory  and  experiment  for  the  first  set  of  measurements  made  in  pure 
oxygen  freestream  with  a  total  enthalpy  of  4  MJ/kg  are  shown  in  Figure  7.  As  with  air,  excellent 
agreement  is  achieved  between  the  heat  transfer  and  pressure  measurements  and  the  predictions 
for  flows  where  the  nonequilibrium  effects  are  small.  However,  when  we  increase  the  total 
enthalpy  level  in  the  freestream  for  a  pure  oxygen  flow  in  Figure  8,  again  the  length  of  the 
separated  region  is  significantly  underpredicted. 
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Figure  7. 


Comparison  of  heat  transfer  and  pressure  measurements  with  Nompelis  Navier-Stokes 
predictions  for  100%  oxygen  freestream  at  a  total  enthalpy  of  4  MJ/kg 
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Figure  8.  Comparison  of  heat  transfer  and  pressure  measurements  with  Nompelis  Navier-Stokes 
predictions  for  100%  oxygen  freestream  at  a  total  enthalpy  of  10  MJ/kg 

Of  significant  interest  in  the  10  MJ/kg  test  case  is  that  the  initial  cone  pressure  significantly  over¬ 
predicted  where  the  flow  is  attached  (as  illustrated  Figure  9(b))  which  suggests  that  the 
calculation  of  the  nozzle  flow  could  be  inadequate.  This  is  a  significant  finding  because  the  same 
models  describing  oxygen  chemistry  and  thermodynamics  that  are  used  on  the  double  cone  are 
also  used  to  predict  the  state  of  the  freestream  after  expanding  through  the  nozzle. 

Decreasing  the  oxygen  content  improves  the  agreement  in  cone  pressure  as  shown  in  Figure  9(c); 
however,  the  heating  rate  on  the  cone  is  over-predicted  for  both  catalytic  and  non-catalytic 
boundary  conditions.  Increasing  the  density  of  the  freestream  flow,  which  implies  an  increase  in 
the  number  of  molecular  collisions  occurring  in  the  expanding  nozzle  flow,  significantly 
improves  the  agreement  as  shown  in  Figure  9(d).  Clearly  more  accurate  models  are  required  to 
describe  the  chemistry  in  the  freestream  before  we  can  evaluate  the  accuracy  of  the  models  of 
flow  chemistry  needed  to  describe  the  shock  layer  flows  over  the  double  cone.  Further  detailed 
flowfield  measurements  are  required  to  resolve  this  problem  and  our  major  emphasis  in  the 
future  will  be  investigating  these  flows  in  the  expansion  tunnel  where  the  conditions  will  be 
selected  to  obtain  flows  with  negligible  nonequilibrium  in  the  freestream. 
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(a)  100%  02/4  MJ/kg 


mimiwuh 


(b)  100%  02/ 10  MJ/kg 


(c)  50%  02  -  50%  Ar  /  10  MJ/kg 


(d)  50%  02  -  50%  Ar  /  10  MJ/kg 


Figure  9.  Variation  of  pressure  and  heat  transfer  on  forecone  with  total  enthalpy  and  freestream 
static  pressure 
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Table  1.  Freestream  conditions  for  double  cone  study 


Run# 

total 

enthalpy 

pressure 

density 

temperature 

vibrational 

temperature 

velocity 

frozen  Mach 

cN2 

cN 

c02 

cO 

cAr 

(MJ/kg) 

(Pa) 

(kg/m3) 

(K) 

(K) 

(m/s) 

78 

3.48 

38.5 

0.0009894 

131 

2,425 

2,456 

10.5 

1.000000 

0.000000 

79 

5.29 

63.4 

1.316E-03 

162 

2,690 

3,073 

11.8 

0.999885 

0.000115 

80 

5.28 

63.6 

1.291E-03 

166 

2,711 

3,067 

11.7 

0.999875 

0.000125 

81 

6.39 

61.3 

9.234E-04 

276 

760 

3,400 

10.8 

0.474165 

0.025835 

0.500000 

82 

8.1 

100.5 

1.045E-03 

391 

660 

3,769 

9.9 

0.454586 

0.045414 

0.500000 

83 

5.12 

35.9 

7.396E-04 

163 

2,849 

3,003 

11.5 

0.999824 

0.000176 

84 

9.27 

109.3 

9.621  E-04 

452 

647 

3,973 

9.6 

0.433079 

0.066921 

0.500000 

85 

10.17 

145.7 

9.243E-04 

586 

687 

4,088 

8.5 

0.665478 

0.084522 

0.250000 

87 

9.85 

164.8 

9.422E-04 

626 

712 

4,019 

8.1 

0.924517 

0.075483 

88 

8.78 

165.3 

1.061E-03 

570 

698 

3,853 

8.2 

0.948209 

0.051791 

90 

3.99 

90.7 

1.834E-03 

190 

1,001 

2,731 

10.4 

0.998556 

0.001444 

91 

10.26 

311.5 

1.550E-03 

729 

773 

4,148 

7.8 

0.938891 

0.061109 

92 

8.77 

205.0 

1.722E-03 

489 

663 

3,960 

9.4 

0.463152 

0.036848 

0.500000 

High-Fidelity  Simulations  of  Turbulent  Flows:  DNS,  WM-LES  and  DES 

A  significant  new  capability  for  the  simulation  of  turbulent  flows  was  developed  under  the 
AFOSR  Center  of  Excellence  support.  Conventional  turbulent  flow  simulation  approaches  use 
Reynolds-averaged  Navier-Stokes  (RANS)  methods,  in  which  the  unclosed  turbulent  transport 
terms  are  modeled  using  a  turbulent  viscosity  and  other  effective  turbulent  transport  coefficients. 
The  models  are  designed  to  work  for  simple  canonical  flows  and  are  then  applied  to  complex 
practical  problems.  The  RANS  approach  is  known  to  work  well  for  certain  classes  of  problems 
such  as  attached  subsonic  boundary  layer  flows.  However,  it  is  seriously  deficient  for  separated 
flows  and  recent  experimental  data  obtained  under  Center  support  at  CUBRC  shows  that  two 
widely-used  RANS  models  do  not  produce  the  correct  heat  flux  for  turbulent  hypersonic  flows. 


As  discussed  above,  a  new  numerical  method  was  developed  that  significantly  reduces  numerical 
dissipation  for  compressible  flows.  This  method  enables  a  new  class  of  simulations  of  turbulent 
high-speed  flows.  Namely,  because  the  method  does  not  artificially  damp  out  small-scale 
motion,  it  can  be  used  for  simulations  that  resolve  many  (if  not  all)  of  the  time-scales  and  length- 
scales  of  turbulent  flows.  This  approach  has  been  applied  to  base  flows  of  re-entry  vehicles,  and 
to  fuel-air  mixing  represented  by  a  jet  or  jets  in  a  supersonic  cross-flow.  This  approach  was  also 
used  to  perform  simulations  of  three  non-trivial  flows  undergoing  transition  to  turbulence  [81, 
88,  102,  104].  Ongoing  work  extends  the  approach  to  the  direct  numerical  simulations  of 
turbulent  boundary  layers,  in  which  all  scales  are  resolved  by  the  simulation. 
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In  the  jet  in  cross-flow  simulations,  the  approach  used  is  a  so-called  wall-modeled  large-eddy 
simulation  or  WM-LES.  Here  a  RANS  model  (in  this  case  the  compressibility-corrected  version 
of  the  Spalart-Allmaras  model)  is  used  near  solid  surfaces  to  provide  a  boundary  condition  to  a 
large-eddy  simulation  (LES)  that  represents  the  large-scale  turbulent  motion  away  from  the 
surface.  This  approach  has  been  shown  to  be  very  effective  for  certain  types  of  flows,  and  we 
have  proven  that  it  can  reproduce  the  available  experimental  data  for  jets  in  supersonic  cross- 
flow  [21,  25,  28,  38,  55,  57,  71,  74,  89,  96,  99,  101,  108].  This  is  very  important  result  because  it 


paves  the  way  for  high-fidelity  simulations  of  full  combustors;  work  is  now  underway  to  use  this 
approach  for  a  representative  geometry  with  multiple  injectors  at  flight-relevant  conditions. 
Figure  10  shows  an  image  from  a  recent  WM-LES  of  a  rectangular  duct  with  5  angled  injectors 
operating. 


Figure  10.  Hydrogen  injection  into  a  supersonic  flow  in  a  rectangular  duct  with  5  angled  (30°) 
injectors  (mole  fraction  of  H2  in  color,  greyscale  represents  density  gradient  magnitude). 
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An  additional  capability  that  the  low-dissipation  numerical  method  enables  is  related  to  transition 
to  turbulence  in  hypersonic  flows.  As  discussed  above,  the  STABL  code  and  other  stability 
analysis  tools  can  correlate  data  for  certain  classes  of  problems.  But  there  are  many  problems  in 
transition  physics  that  are  not  suitable  to  conventional  stability  analyses.  Under  Center  support, 
we  have  applied  the  low-dissipation  approach  to  the  direct  simulation  of  transition  to  turbulence 
for  two  problems:  acoustic  wave  interactions  with  the  Mach  8  flow  over  an  elliptic  cone  at  zero 
angle  of  attack  and  the  instability  growth  in  the  wake  of  a  discrete  roughness  element  in  a  Mach 
6  laminar  boundary  layer.  These  problems  would  not  have  been  feasible  without  the  low- 
dissipation  numerical  method  and  the  work  done  under  Center  support  to  allow  the  code  to  scale 
to  large  numbers  of  processors  and  grid  elements.  Typical  simulations  use  1000  to  2000  cores 
and  require  270  million  element  grids  [88,  102]. 


Figure  1 1 .  Isosurfaces  of  Q  colored  by  velocity  gradient  magnitude  in  the  wake  of  a  cylindrical 
roughness  element  in  a  Mach  6  flow. 

Figure  1 1  shows  a  result  from  the  discrete  roughness  element  simulation.  Flere,  the  roughness 
element  is  located  at  the  lower  left  of  the  image,  and  is  a  cylindrical  geometry  protruding  from 
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the  curved  surface  of  a  wind  tunnel  (the  Purdue  Mach  6  Quiet  Tunnel).  The  variable  plotted  is  Q, 
which  is  the  second  invariant  of  the  velocity  gradient  tensor,  and  represents  the  location  of 
vortical  motion.  The  isosurfaces  of  Q  are  colored  by  the  local  velocity  magnitude.  Note  that 
several  vortices  wrap  around  the  roughness  element  and  are  then  swept  downstream.  These 
interact  with  one  another  and  ultimately  cause  the  wake  to  break  down  to  turbulent  flow. 
Analysis  and  comparison  with  experimental  measurements  made  at  Purdue  show  that  the  wake 
has  a  dominant  unsteadiness  frequency  (about  21  kHz  in  this  case).  We  found  that  this  dominant 
frequency  is  driven  by  the  unsteady  motion  in  the  separation  region  ahead  of  the  protuberance. 

Advanced  Particle-Based  Simulation  Methods  for  Hypersonic  Rarefied  Flows 

AFOSR  Hypersonics  Center  funding  was  used  in  Prof.  Schwartzentruber’s  group  for  continued 
code  development  of  a  state-of-the-art  particle  simulation  method  for  non-equilibrium  gas  flows. 
The  method  is  called  direct  simulation  Monte  Carlo  (DSMC)  and  some  sample  results  are  shown 
below.  DSMC  simulates  the  Boltzmann  equation  and  is  therefore  accurate  for  flows  ranging 
from  continuum  to  free-molecular.  The  method  is  used  to  predict  drag  and  heating  environments 
for  high-speed,  high-altitude  objects  (satellites,  space  vehicles,  etc.);  regimes  where  the 
continuum  fluid  equations  become  inaccurate.  DSMC  moves  millions  of  representative 
molecules  through  a  computational  grid  and  performs  both  gas-phase  collisions  and  gas-surface 
collisions,  including  chemical  reactions.  Due  to  rapidly  increasing  computer  resources  the 
method  is  now  tractable  for  full-scale  engineering  problems,  and  each  year  its  capabilities 
increase.  DSMC  holds  great  potential  for  molecular-based  physics  for  fluid  dynamic  modeling 
involving  chemical  reactions  and  gas-surface  interactions. 

Specifically,  funding  was  used  to  develop  both  adaptive  mesh  refinement  (AMR)  and  ‘cut-cell’ 
capability  for  complex  geometry  into  the  DSMC  code.  The  code  has  been  threaded  using 
OpenMP  for  parallel  simulations  on  shared-memory  architectures.  Finally,  a  preliminary  MPI 
version  has  been  implemented  for  parallel  simulations  on  distributed  memory  architectures  [75, 
76,  92,  93,  98], 

In  order  to  obtain  accurate  CFD  solutions  for  hypersonic  flows,  grids  must  be  aligned  with  strong 
gradients  in  the  flow,  often  requiring  an  unstructured  grid  and  a  user  with  advanced  grid- 
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generation  skills.  In  contrast,  accurate  DSMC  solutions  can  be  obtained  on  Cartesian  grids.  This 
offers  a  number  of  advantages  for  large  parallel  simulations  however,  requires  “cut-cell” 
algorithms  to  cut  complex  geometry  surfaces  from  a  Cartesian  flow  field  grid.  Finally,  for 
accuracy  and  efficiency,  DSMC  grid  cells  must  be  adapted  to  the  local  mean  free  path 
(proportional  to  the  local  density).  Essentially  the  optimal  DSMC  grid  is  coupled  to  the  final 
solution  itself  and  thus  naturally  requires  AMR.  The  key  idea  is  that  both  AMR  and  cut-cell 
procedures  can  be  fully  automated  and  parallelized  in  DSMC  (Cartesian  grids),  potentially 
eliminating  the  user-time  for  grid  generation.  The  DSMC  code  at  Minnesota  is  now  able  to 
automatically  generate  and  refine  a  three-level  Cartesian  flow  field  grid  around  complex  3D 
geometries  and  obtain  accurate  solutions  with  no  user  intervention.  An  example  solution  for  a 
satellite  geometry  is  shown  below. 


Figure  12.  3D  DSMC  solution  for  hypersonic  flow  over  a  satellite  in  low  Earth  orbit, 
demonstrating  complex  geometry  capability. 


With  the  mechanics  of  the  DSMC  code  complete  (i.e.  surface  geometry,  flow  field  grid,  particle 
movement  and  sorting  algorithms),  the  DSMC  method  becomes  highly  modular  with  respect  to 
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physical  modeling.  Specifically,  during  a  single  time  step,  all  chemistry  models  operate  only  on 
the  particles  within  each  cell  and  do  not  depend  on  neighboring  cells.  This  enables  rapid 
incorporation  of  new  physical  models  into  the  large-scale  Minnesota  DSMC  code.  Center 
funding  has  partially  supported  ab-initio  computational  chemistry  research  on  gas-surface  and 
gas-phase  collisions.  Specifically,  Prof.  Schwartzentruber’s  group  has  developed  a  novel 
accelerated  Molecular  Dynamics  (MD)  method  specifically  for  dilute  gases.  This  method  enables 
the  efficient  simulation  of  normal  shock  waves  using  only  atoms  and  associated  inter-atomic 
potentials  (no  collision  models,  transport  closures,  or  rate  models).  The  method  has  been 
validated  with  experimental  data  for  shock  waves  in  argon  as  well  as  argon-xenon  and  argon- 
helium  mixtures.  The  data  produced  by  such  first-principles  MD  simulations  is  then  used  to 
improve  DSMC  collision  models  used  in  the  large-scale  Minnesota  code.  Ongoing  work  focuses 
on  polyatomic  gases  including  rotation/vibration  excitation  and  dissociation  of  air  in  shock  layers 
[64,  65,  67], 
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Additional  Specific  Topics  Addressed  Under  Center  Support 

In  the  following,  additional  details  concerning  specific  sub-projects  accomplished  under  Center 
support  are  listed  in  approximate  chronological  order.  Most  of  this  work  is  supported  by  detailed 
publications  listed  in  this  report. 


1.  A  new  Mach  10  scramjet  engine  was  designed  by  Center  personnel  and  successfully 
tested  in  the  CUBRC  LENS-I  facility  in  support  of  the  DARPA  HyCAUSE  program. 
This  engine  uses  an  inward-turning  inlet  design  and  is  the  first  such  engine  ever  tested  at 
full  scale  a  fully  duplicated  flight  conditions.  Center  personnel  and  computational  tools 
developed  at  the  Center  were  instrumental  in  the  robust  performance  obtained  from  this 
engine.  This  inlet/engine  concept  was  flight-tested  at  the  Australian  range  under  support 
from  the  HyCAUSE  program;  unfortunately  the  flight  test  was  largely  unsuccessful 
because  of  problems  with  orienting  the  flight-test  payload  during  re-entry. 

2.  The  AFOSR  STAR  program,  which  was  partly  funded  by  Center  support,  made  great 
progress  in  transitioning  high-fidelity  mechanism-based  stability  analysis  tools  to  the 
wider  community.  The  STABL  tool  (Stability  and  Transition  Analysis  for  hypersonic 
Boundary  Layers)  has  undergone  significant  improvements  and  is  now  in  use  at  many 
academic  institutions  and  government  labs.  STABL  was  used  to  predict  the  transition 
location  on  the  Lockheed-Martin  STAR- 1 A  geometry  to  support  the  development  of  the 
DARPA  FALCON  program  HTV-1  vehicle.  We  supported  the  FALCON  HTV-2 
aero  thermal  and  transition  analysis,  and  work  is  underway  to  analyze  transition  on  the  X- 
51  A.  STABL  is  being  used  to  support  the  Raytheon  KEI  development  program.  STABL 
was  also  used  to  study  transition  to  turbulence  on  the  HIFiRE-1  and  -5  flight  vehicles  to 
design  the  experiments  [12,  24,  34,  46-49,  50,  69]. 

3.  Computational  tools  developed  at  the  Center  were  used  to  support  the  design  of  the 
DARPA  FALCON  Hypersonic  Cruise  Vehicle  (HCV)  inlet  and  engine  flow  path. 

4.  A  new  low  dissipation  method  for  the  large-eddy  simulation  of  compressible  turbulence 
has  been  developed.  This  approach,  which  explicitly  conserves  kinetic  energy,  captures  a 
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much  larger  range  of  length  scales  than  conventional  upwind  methods.  It  will  enable  the 
LES  of  compressible  turbulent  mixing  and  combustion  [63]. 

5.  A  novel  shock-capturing  scheme  was  developed  for  unstructured  grids;  the  predictor- 
corrector  approach  is  applied  in  the  vicinity  of  detected  shock  waves.  It  has  the  potential 
to  improve  solution  accuracy  and  robustness  for  high-Mach  number  flows,  particularly 
for  direct  numerical  simulations  (DNS)  and  large-eddy  simulations  (LES). 

6.  A  new  parallel  implicit  unstructured  grid  code  was  developed  and  validated  [5,  1 1].  This 
code  uses  the  efficient  data-parallel  line  relaxation  approach  for  unstructured  grids, 
combining  the  efficiency  of  structured  grid  approach  with  the  flexibility  of  unstructured 
grids  for  complex  geometry  applications.  This  code,  US3D,  is  becoming  widely  used;  it 
is  now  in  use  at  Kirtland  AFB  in  the  AFRL  Space  Vehicles  Directorate  and  at  Wright- 
Patterson  AFB  in  the  AFRL  Air  Vehicles  Directorate,  as  well  as  at  NASA  and  several 
academic  institutions  and  industry  partners. 

7.  A  nozzle  flow  analysis  tool  was  developed  and  released  through  the  AEDC  [14].  This 
code  was  used  to  analyze  the  AEDC  Tunnel  9  flow  and  recent  double-cone  experiments 
were  simulated  at  Tunnel  9  flow  conditions.  Center  researchers  collaborated  closely  with 
AEDC  Hypervelocity  Tunnel  9  personnel  to  simulate  the  Mach  14  nozzle  with  the  goal 
of  understanding  its  non-ideal  performance.  As  part  of  this  work  a  novel  approach  for 
generating  grids  in  hypervelocity  flows  was  conceived,  in  which  the  grid  lines  are 
aligned  with  direction  of  the  characteristics  in  the  supersonic  portion  of  the  flow.  This 
approach  was  shown  to  resolve  small  disturbances  much  more  accurately  then 
conventional  grids  that  are  not  aligned  with  those  weak  disturbances.  This  approach  has 
the  potential  to  greatly  improve  the  accuracy  of  wind  tunnels  designs,  and  more  broadly, 
the  simulation  of  supersonic  flows  in  general  [107]. 

8.  New  code  validation  experiments  were  conducted  and  the  new  non-intrusive  diagnostics 
for  nitric  oxide  and  water  vapor  concentration  were  used  to  help  quantify  the  flow  state 
in  the  LENS-I  test  section  [22,  40].  These  are  the  first  nitric  oxide  measurement  ever 
made  in  the  test  section  of  a  hypervelocity  facility.  Work  was  conducted  with  the  attempt 
to  explain  the  very  puzzling  results  obtained  from  the  non-intrusive  diagnostics 
measurements.  This  work  has  resulted  in  a  series  of  additional  papers;  unfortunately,  this 
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problem  remains  unsolved  and  will  be  the  subject  of  part  of  the  proposed  computational 
and  experimental  program.  The  present  thinking  is  that  the  modeling  of  oxygen 
recombination  in  high-enthalpy  facility  nozzles  is  significantly  inaccurate,  which  leads  to 
the  discrepancies  between  simulation  and  experiment.  It  will  take  considerably  more 
effort  to  quantify  this  effect  and  to  develop  predictive  models  of  the  recombination 
process.  Furthermore,  analysis  shows  that  the  electronically  excited  states  may  play  a 
role  in  the  recombination  process,  which  will  require  a  complete  refonnulation  of  all 
modeling  approaches  for  high-temperature  flows  [80,  85,  97]. 

9.  Extensive  experimental  programs  were  conducted  in  the  CUBRC  facilities,  including 
studies  on:  high-enthalpy  flows  over  spheres,  cylinders,  capsules,  and  double-cones; 
catalytic  heating  on  capsule  configurations;  tripping  and  transition  on  scramjet  engine 
inlets.  Experiments  were  performed  to  support  the  AFRL  HiFIRE  Flight  1.  Extensive 
simulations  were  also  perfonned  to  support  these  flight  experiments.  In  addition,  studies 
were  made  to  produce  comparisons  of  predicted  transition  locations  on  the  HIFiRE-1 
geometry  with  measurements  in  the  NASA  Langley  Mach  6  facility.  The  key  results  of 
these  studies  are  that  the  eN  approach  correlates  the  HIFiRE-1  configuration  data  across 
two  very  different  facilities  and  very  different  flow  field  conditions  34,  46-48,  50].  This 
is  an  important  result  and  shows  the  potential  of  the  mechanism-based  stability  and 
transition  analysis  approach. 

10.  Under  Center  support  there  were  modifications  and  improvements  to  the  LENS-X 
expansion  tunnel  for  real  gas  studies  at  velocities  of  10,000  to  27,000  ft/s  with 
significantly  reduced  issues  of  test  gas  reaction.  Studies  were  conducted  to  make  cross¬ 
facility  comparisons;  this  work  demonstrated  that  the  expansion  tunnel  is  a  viable  facility 
for  producing  extremely  high-enthalpy  flows  with  very  low  freestream  contamination 
and  thermo-chemical  nonequilibrium.  Subsequent  studies  have  confirmed  this  result  with 
the  new  CUBRC  LENS-XX  expansion  tunnel  facility. 

1 1.  A  Workshop  on  Hypersonics  was  held  in  Minneapolis  September  14-15,  2005. 

12.  There  were  extensive  interactions  with  industry,  including  AeroSystems  Engineering, 
Alliant  Tech  Systems,  Honeywell,  Boeing,  Lockheed-Martin,  Pyrodyne,  Aerojet, 
Northrop-Grumman,  Raytheon  Missiles  and  Space,  and  GoHypersonic  Inc. 
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13.  Many  papers  have  been  published  (approximately  100  conference  papers  and  35  journal 
articles);  numerous  invited  and  contributed  talks  were  given  by  Center  researchers. 
Approximately  fifteen  Ph.D.  and  30  M.S.  students  graduated  from  the  University  of 
Minnesota. 

14.  The  University  of  Minnesota  course  in  Hypersonic  Aerodynamics  was  extensively 
updated.  The  course  notes  are  being  put  into  the  format  for  a  textbook  on  the  subject.  A 
new  course  on  Molecular  Gas  Dynamics  was  developed  by  Prof.  Schwartzentruber  to 
teach  the  key  aspects  of  the  Direct  Simulation  Monte  Carlo  (DSMC)  method  and 
concepts  in  rarefied  gas  dynamics  with  application  to  hypersonic  flows  at  high  altitude 
conditions. 

15.  Center  personnel  have  served  on  several  advisory  committees  and  panels  for  NASA, 
DoD,  industry,  Sandia  National  Laboratories,  and  the  Air  Force. 

16.  The  Center  has  been  providing  support  for  the  HIFiRE  Flight  5  vehicle  design;  this  has 
involved  close  collaborations  with  AFRL  personnel  and  resulted  in  publications 
documenting  the  work. 

17.  Experimental  studies  have  been  conducted  in  the  LENS  II  tunnel  to  examine  the 
aerothermal  characteristics  of  a  scramjet  powered  hypersonic  vehicle  designed  for  Mach 
6  flight.  Measurements  made  in  these  studies  are  being  employed  to  examine  and 
improve  the  turbulence  modeling  used  in  the  DPLR  and  US3D  numerical  prediction 
codes,  as  well  as  other  standard  CFD  codes. 

18.  Extensive  experimental  studies  were  conducted  in  the  LENS  tunnels  in  support  of  the 
AFOSR  HIFiRE  and  the  Air  Force  FALCON  program.  Measurements  made  in  these 
programs  were  employed  to  examine  the  prediction  of  transition  on  hypersonic  vehicles 
and  to  develop  improved  models  of  turbulence  for  predicting  regions  of  shock 
wave/boundary  layer  interaction. 

19.  Experimental  studies  have  continued  in  the  LENS  I  hypervelocity  tunnel  to  examine  the 
nonequilibrium  characteristics  of  a  high-enthalpy  flowfield.  Analysis  of  these 
measurements  is  continuing  with  the  support  of  results  from  numerical  prediction 
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techniques  including  sophisticated  models  of  flowfield  chemistry.  Much  of  this  work 
focused  on  the  double  cone  configuration. 

20.  Significant  efforts  were  made  in  the  design  and  evaluation  of  the  new  LENS  XX 
expansion  tunnel  designed  to  generate  velocities  from  10,000  to  27,000  ft/s.  In  concert 
with  the  design  of  this  facility,  detailed  simulations  were  made  to  optimize  a  contoured 
nozzle  which  will  meet  the  requirements  of  the  LENS  XX  expansion  tunnel.  The  LENS 
XX  facility  represents  a  completely  new  ground-test  capability  for  high-enthalpy  flows. 
Key  aspects  are  that  many  re-entry  flight  conditions  can  be  replicated  in  the  facility,  and 
that  the  expansion  tunnel  has  very  low  levels  of  chemical  reaction  and  internal  energy 
excitation  are  present  in  the  freestream  test  gas.  This  is  crucial  for  high-fidelity  testing, 
radiative  emission  measurements,  and  validation-quality  data  for  thermo-chemical  model 
validation  and  calibration. 

21.  The  Center’s  computational  tools  were  used  to  support  the  design  of  the  DARPA  ASET 
inward-turning  inlet.  The  inlet  was  optimized  and  studied  under  a  wide  range  of  flow 
conditions  using  Center  developed  simulation  tools;  this  inlet  and  flow  path  was  tested  at 
the  AEDC  Tunnel  9  facility  and  potentially  has  application  to  future  hypersonic  cruise 
missile  concepts. 

22.  The  University  of  Minnesota  computational  fluid  dynamics  code,  US3D,  was  extended 
to  larger  and  more  complex  problems;  a  recent  example  simulation  used  128  million 
hexahedral  elements  and  a  5 -species  finite-rate  thermo-chemical  model  for  the  Space 
Shuttle  Orbiter  flow  field.  To  the  best  of  our  knowledge,  this  is  the  largest 
aerothermodynamics  simulation  ever  performed.  Other  large-scale  simulations  are  being 
done  to  support  AF  programs  (HIFiRE-5  and  X-51A,  for  example). 

23.  A  possible  explanation  for  the  observed  differences  in  reflected  shock  tunnel 
perfonnance  and  simulations  was  proposed;  this  involves  a  new  mechanism  in  which 
electronically  exited  states  of  oxygen  and  carbon  dioxide  are  formed  during  the 
recombination  process  in  the  facility  nozzle.  New  work  is  progressing  to  understand  and 
quantify  this  effect  [80,  85,  92,  97]. 

24.  The  US3D  computational  fluid  dynamics  code  was  extended  to  include  several  new  time 
integration  methods  required  for  high  order  of  accuracy  simulations.  Several  new 
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approaches  were  investigated  and  one  novel  method  was  developed  that  shows  great 
promise  for  wall-modeled  large-eddy  simulations  (WM-LES).  This  approach  makes 
minor,  but  critical,  modifications  to  the  classic  second-order  Crank-Nicholson  method 
that  provides  stability  and  accuracy  for  unsteady  hypersonic  flows.  This  method  has  been 
shown  to  permit  large  implicit  time  steps  in  the  near-wall  RANS  region,  while 
maintaining  accuracy  and  low  dissipation  levels  in  the  LES  region.  This  approach  is 
much  more  cost-effective  than  competing  methods,  at  least  for  the  certain  classes  of 
flows  studied  under  Center  support.  For  example,  recent  simulations  of  jets  in  supersonic 
cross-flow  are  enable  by  this  approach. 

25.  The  basic  implicit  time  integration  method  in  US3D  was  extended  from  a  data-parallel 
line-relaxation  (DPLR)  approach  to  a  Generalized  Method  of  Residual  (GMRes)  method 
to  see  if  it  is  possible  to  improve  convergence  rates  on  non-ideal  grids  and  for  grids  that 
are  not  suited  to  the  line-relaxation  method.  It  was  found  that  under  some  conditions  the 
GMRes  method  is  competitive  with  the  DPLR  method,  however  this  very  much  depends 
on  the  type  of  pre-conditioner  used,  the  optimization  of  settings  for  the  GMRes  method 
and  other  details.  In  general,  the  DPLR  approach  works  well  and  in  some  cases  can  be 
made  to  be  significantly  faster  than  GMRes  by  adjusting  the  number  of  sub-iterations 
used  during  the  relaxation  step  [44,  103]. 

26.  Related  to  the  GMRes  approach,  a  Newton-like  time  integration  method  was  also  studied 
with  the  goal  of  reducing  costs  of  performing  shape  optimization  on  hypersonic  inlets, 
for  example.  It  was  found  that  by  converting  the  DPLR  method  to  a  Newton  method  and 
using  a  large  number  of  sub-iteration  steps,  it  is  possible  to  reduce  convergence  times  by 
a  substantial  margin  for  cases  in  which  the  flow  conditions  or  geometry  do  not  change 
much  [94,  100].  For  example,  during  a  shape  optimization,  the  derivatives  of  the 
objective  function  with  respect  to  geometry  variables  are  required.  This  requires  a 
baseline  CFD  simulation  and  then  a  series  of  simulations  on  slightly  perturbed  grids. 
This  Newton  approach  can  significantly  reduce  the  cost  of  the  secondary  simulations  (by 
approximately  a  factor  of  10),  making  shape  optimization  much  less  computationally 
demanding.  Furthermore,  this  approach  can  be  used  to  rapidly  investigate  aerodynamic 
perfonnance  as  a  function  of  flight  conditions  and/or  angle  of  attack  by  propagating  a 
solution  forward  from  a  baseline  simulation. 
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27.  Significant  and  long-term  support  was  provided  through  the  Center  for  experiments  at 
the  Caltech  T5  Free-Piston  Shock  Tunnel.  Center  researchers  performed  numerous 
computations  and  interacted  closely  with  Caltech  and  AFRL  personnel  to  understand  the 
experimental  data  and  to  design  experiments  [61,  82].  The  goal  of  this  work  was  to  inject 
carbon  dioxide  into  a  hypervelocity  flow  to  control  transition  to  turbulence.  The  Center 
simulated  the  experiments,  performed  many  STABL  analyses  of  the  boundary  layers  to 
understand  how  the  CO2  affected  the  flow  stability  and  how  the  injection  process  could 
interfere  with  transition  delay  by  destabilizing  the  boundary  layer.  The  Center  work  was 
extremely  important  in  helping  characterize  the  flow  physics  and  ultimately  helped  lead 
to  the  success  of  the  experiments. 

28.  Close  collaborations  between  Center  researchers  and  Purdue  University  personnel 
resulted  in  a  new  understanding  of  wake  transition  physics  due  to  discrete  roughness 
elements  and  protuberances.  In  particular,  detailed,  highly-resolved  simulations  of  the 
Purdue  Mach  6  Quiet  Tunnel  experimental  configuration  were  carried  out  by  Center 
researchers  [88,  102].  Comparisons  with  experimental  data  were  made;  the  computations 
indicated  that  the  source  of  the  dominant  frequency  of  wake  unsteadiness  originates  in 
the  region  upstream  of  the  discrete  roughness  element.  Experimental  pressure 
transducers  were  then  placed  in  that  region  and  were  found  to  measure  very  similar 
disturbance  properties  to  that  predicted  by  the  simulations.  This  is  probably  the  first  time 
that  a  direct  numerical  simulation  was  used  to  make  a  prediction  about  a  realistic  and 
technically  relevant  high-speed  flow  field.  This  work  would  not  have  been  possible 
without  the  numerical  methods  and  code  development  performed  under  Center  support. 

29.  Similar  to  the  above,  the  US3D  code  was  used  to  study  cross-flow  transition  on  a  cone  at 
conditions  relevant  to  the  Purdue  Mach  6  Quiet  Tunnel  flow  [105],  Simulations 
produced  cross-flow  disturbances  that  are  similar  to  those  observed  in  wind  tunnel 
measurements.  This  work  is  still  preliminary,  but  shows  the  potential  for  direct 
simulations  of  this  type  of  transition  mechanism  in  hypersonic  flows. 

30.  A  very  large  body  of  work  was  conducted  related  to  wall-modeled  large-eddy 
simulations  (WM-LES),  with  emphasis  on  fuel-air  mixing  simulations  at  technologically 
relevant  conditions  [21,  25,  28,  30,  38,  55,  57,  71,  74,  89,  96,  99,  101,  108].  This  work 
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has  shown  that  with  the  use  of  appropriate  low-dissipation  high-order  of  accuracy 
numerical  methods,  large  high-quality  carefully  designed  grids,  and  certain  variants  of 
turbulence  models  and  inflow  conditions,  it  is  possible  to  reproduce  experimental  data 
for  fuel-air  mixing.  This  is  an  extremely  important  conclusion  because  it  paves  the  way 
for  detailed  simulations  of  actual  combustor  configurations,  and  ultimately  for  the 
optimization  of  hypersonic  air-breathing  propulsion  system  flow  paths.  This  is  a  major 
accomplishment  of  the  Center  and  has  the  potential  to  revolutionize  the  efficiency  and 
robustness  of  future  propulsion  systems.  We  are  now  working  to  transition  this  capability 
to  the  AFRL  and  to  potential  industry  users. 

3 1 .  Related  to  the  above  accomplishment,  Center  funds  were  used  to  develop  an 
experimental  configuration  for  obtaining  high-quality  open-literature  data  for  code 
validation  of  fuel-air  mixing  and  combustion  at  relevant  conditions  and  with  multiple 
interacting  injection  ports.  This  “combustion  duct”  was  designed  to  have  some  similarity 
to  the  AFRL  HIFiRE-2  flight  experiment  configuration,  but  is  completely  generic  and  is 
not  in  any  way  optimized  to  produce  a  high-quality  flight  system  relevant  flow.  It  uses 
simple  angled  round  injection  ports  and  can  be  configured  in  many  different  ways  to 
allow  specific  data  to  be  obtained.  Up  to  9  individual  fuel  injection  ports  may  be 
operated.  The  experimental  rig  has  optical  access  that  will  enable  laser-based  diagnostics 
to  be  used  to  characterize  the  flow  field.  Initial  experiments  were  run  at  CUBRC  and 
new  experiments  are  planned  with  additional  diagnostics  methods  in  place  to  extend  the 
value  of  the  data.  At  the  University  of  Minnesota,  the  CUBRC  combustion  duct  has  been 
simulated  using  the  WM-LES  approach  with  various  injector  configurations  in  place, 
and  at  several  flow  conditions.  Initial  comparisons  with  the  limited  experimental  data  are 
very  encouraging  [101].  If  it  is  possible  to  obtain  agreement  between  simulation  and 
experiment  for  this  configuration,  it  will  provide  confidence  that  the  numerical  methods 
and  turbulence  modeling  approaches  can  be  used  for  system-relevant  air-breathing 
propulsion  systems. 

32.  As  the  WM-LES  approaches  and  associated  numerical  methods  have  gained  maturity, 
Center  researchers  have  begun  to  perform  direct  numerical  simulations  (DNS)  of 
hypersonic  turbulent  flows.  This  approach,  which  resolves  all  physically-relevant  length 
and  time  scales,  promises  to  provide  statistically  exact  data  of  important  turbulent  flows. 
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For  example,  a  puzzling  result  has  been  observed  on  the  HIFiRE-1  blunt  cone 
configuration.  Here,  the  measured  turbulent  heat  flux  is  best  predicted  with  a  very  simple 
(Baldwin-Lowmax)  turbulence  model,  while  more  physically-relevant  models  over¬ 
predict  the  measured  heat  flux  by  substantial  margins.  This  is  of  course  a  major  concern 
and  has  serious  implications  for  system  design  uncertainty  and  thennal  protection  system 
sizing.  Through  Center  support,  it  is  now  possible  to  perfonn  DNS  of  these  flows  at 
actual  wind  tunnel  conditions;  though  these  simulations  are  time-consuming  (requiring 
about  2  months  on  1000  processor  cores  with  a  300  million  element  grid)  ,  they  have  the 
potential  to  improve  the  turbulence  models  and  explain  why  the  heat  flux  is  so  poorly 
predicted.  This  capability  should  enable  new  simulations  to  help  explain  many 
outstanding  problems  in  hypersonic  flow  physics. 

33.  Related  to  the  WM-LES  and  DNS  capability,  is  the  area  of  chemistry-turbulence 
interactions.  As  described  above,  it  is  now  possible  to  predict  fuel-air  mixing  in  practical 
supersonic  air-breathing  propulsion  configurations.  However,  a  major  issue  remains: 
how  to  model  the  interaction  of  turbulent  motion  with  chemical  reactions  and  the 
combustion  process  in  general.  A  huge  body  of  work  exists  for  the  modeling  of  this 
interaction  at  low-speed  combustion  conditions,  but  there  is  no  accepted  approach  for 
high-speed  (supersonic)  flow  conditions.  Because  of  the  high  speeds,  the  relative  time 
scales  of  fluid  motion  and  chemical  reaction  are  considerably  different  and  thus  the  well- 
established  models  cannot  be  applied  without  modification.  Under  Center  support,  work 
in  this  area  was  initiated,  but  was  not  completed.  The  US3D  code  was  modified  to  allow 
Lagrangian  particle  tracking  so  that  a  Filtered  Mass  Density  Function  (FMDF)  approach 
can  be  used  to  represent  subgrid  scale  chemistry-turbulence  interactions  in  a  large-eddy 
simulation  approach;  the  FMDF  approach  has  been  very  successful  at  low  speeds  and 
should  be  readily  adaptable  to  the  high-speed  flow  environment.  With  further 
development  this  approach  shows  promise  for  providing  a  rigorous  chemistry-turbulence 
interaction  closure. 

34.  A  novel  accelerated  molecular  dynamics  (MD)  approach  was  developed  that  allows  low- 
density  flows  to  be  simulated  at  a  near  first-principles  level.  The  accelerated  MD 
collapses  the  simulation  time  that  occurs  between  collisions  so  that  the  simulations  do 
not  compute  the  particle  interactions  when  they  are  not  interacting  and  are  on  ballistic 
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trajectories.  This  is  an  important  advance  and  will  enable  new  much  more  detailed 
simulations  of  critical  physical  processes  in  high-temperature  flows.  An  important 
application  is  to  shock  waves  so  as  to  better  understand  the  physics  that  occur  within  the 
highly  nonequilibrium  region  associated  with  a  shock  wave  [64,  65,  67]. 
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AIR  FORCE  RELATED  TRANSITIONS 


The  nozzle  flow  simulation  CFD  code  has  been  transitioned  to  the  Arnold  Engineering  and 
Development  Center  (AEDC)  for  analysis  of  the  flow  in  the  Aero-Propulsion  Test  Unit  (APTU) 
and  Tunnel  9. 

Extensive  simulations  of  the  AEDC  Tunnel  9  nozzle  flow  were  conducted  and  double-cone 
experiments  were  designed  for  the  facility.  Computations  were  performed  to  support  the  analysis 
of  the  double  cone  experiments. 

The  STABL  (Stability  and  Transtion  Analysis  for  hypersonic  Boundary  Layers)  code  has  been 
transitioned  to  many  sites  including  AFRL  Air  Vehicles,  Wright-Patterson  AFB.  The  code  is  in 
active  use  at  these  institutions  for  hypersonic  transition  analysis  and  comparison  to  flight  and 
ground  test  data,  including  several  contractors  that  are  using  STABL  for  AF  and  other  DoD 
work. 

The  US3D  code  has  been  transitioned  to  AFRL  Space  Vehicles,  Kirtland  AFB  and  AFRL  Air 
Vehicles,  Wright-Patterson  AFB. 

There  were  extensive  interactions  and  collaborations  between  Center  researchers  and  Caltech 
and  AFRL  Propulsion  Directorate,  Edwards  AFB. 

Graduate  students  and  researchers  spent  time  during  the  summer  at  Wright-Patterson  AFB 
working  with  AFRL  personnel. 
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